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Abstract 



Motivated by numerical experiments and studies of quantum systems 
which are classically chaotic, we take a Random Matrix description of a Hard- 
sphere gas to Statistical Mechanical description. We apply this to Anyon gas 
and obtain a formal expression for the momentum distribution. Various lim- 
iting situations are discussed and are found in agreement with the well-known 

results on Hard-sphere gas in low-density regime. 
PACS numbers:05.30.-d 05.45. +b 
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Random Matrix Theory and Statistical Mechanics are employed for the description of 
many-body systems like nuclei, metallic clusters and so on. During the last decade or so, it 
has become increasingly clear that spectral fluctuations of simple quantum systems [[[J whose 
classical counterpart are chaotic, as well as the eigenfunctions, have features |2j in striking 
agreement with random matrix theory (RMT). Originating from a conjecture of Berry |3j 
regarding the eigenfunctions as a random superposition of plane waves, it was shown that 
a contact with statistical mechanics can be realised [[§]. In this Letter, we show that this 
connection is indeed a way to relate RMT and Statistical Mechanics. 

It is well-known that statistical mechanics in two dimensions is at the core of several 
phenomena of great interest such as the fractional Hall effect, high-T c superconductivity, 
and others ]5|||. Unlike in three dimensions the analysis here is plagued with difficulties 
stemming from an unusual, non-local exchange interaction potential. As a result identi- 
cal particles in two dimensions exhibit statistics interpolating between the Fermi-Dirac and 
Bose-Einstein distributions, generically these are termed as Anyons 0. In this Letter, we 
address this fundamental problem in two dimensions and cast it systematically employing 
an ansatz for the eigenfunctions consistent with the connection mentioned above, and in- 
corporate the Braid group governing exchange symmetries. We show a way to arrive at 
the momentum distribution upto 0(h 2 ) and express it in terms of quantities involving the 
counting problem in Braid group. Statistical mechanics of Anyons has been studied in the 
past H and it is very interesting to see that the second virial coefficient obtained then and 
now, give rise to a Sum Rule. However, in marked distinction with previous works, calcula- 
tion of higher-order virial coefficients is also possible here - a facet which makes the present 
approach a novel alternative. To this end, we argue that the third virial coefficient has the 
form in agreement with some recent results. We re-iterate that the aim of this Letter is to 
arrive at quantum statistical mechanics without having to make a hypothesis about thermal 
bath, and that the present approach is guided by random matrix theory. 

The fundamental hypotheses of RMT are (a) the hamiltonian of a given system belongs 
to an ensemble of hamiltonians, and (b) real quantum mechanics is enough for description of 
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physical systems if time-reversal symmetry is taken into account appropriately ||. It is the 
latter of the two that leads to exactly three universality classes in the RMT - Orthogonal 
(even spin, time-reversal invariant, OE) Unitary (time-reversal non-invariant, UE) and Sym- 
plectic Ensembles (odd spin, Time- Reversal Invariant, SE) |Kj. Owing to the spin-statistics 



connection in two dimensions [II], we know that there is fractional statistics and fractional 
spin. In the context of symmetries, in two dimensions, parity (P) and time-reversal (T) both 
are broken. 

Let us consider a system of N hard spheres ('discs' in two dimensions), each of radius a, 
enclosed in a box of edge-length L + 2a. Centres of two hard spheres xi and Xj are such that 
\%i > 2a. The canonical pair of coordinates describing these particles are (X, P) where 
X = (xi,x 2 ,---,Xn), P = (pi,P2, • • ■ ,Pn)- Energy eigenfunctions, i[) a (X) corresponding to 
eigenvalue E a vanish on the boundary of the enclosure. A typical eigenfunction is irregular, 
with a Gaussian amplitude distribution and the spatial correlation function of the same is 
consistent with the conjecture of Berry which allows us to represent this eigenfunction as a 
superposition [[| : 

MX) = N a J d dN PA a (P)5(P 2 - 2mE a ) e** p (1) 

with N a given by the normalization constant, and A' a s satisfying the two-point correlation 
function 

5 dN (P-P') 

d denotes the number of coordinate-space dimensions. The average in @ is a matrix- 
ensemble (ME) average which originates from the fact that the hamiltonian, H of the system 
belongs to an ensemble of matrices satisfying associative division algebra |]I0 , |I2]| in consis- 
tency with quantum mechanics. The eigenstate ensemble (EE) used in [|J is nothing but a 
consequence of underlying matrix ensemble in RMT, the eigenfunctions then satisfy all the 



properties numerically observed and analytically represented in ([!]), (^) []I3"| |. The correla- 
tion functions @ decide whether time-reversal symmetry is preserved {A* a (P) = A a (—P)) 



— * — * 

or broken (A* a (P) ^ A a (—P)), accordingly the corresponding matrix ensemble belongs to 
OE or UE respectively. As noted in the higher-order even-point correlation functions 
factorize and the odd-ones vanish. A very important aspect of the ansatz (fj), (||D is that the 
Wigner function corresponding to ip a (X) is microcanonical, or, is proportional to S(H — E a ) 
which, in a sense, incorporates ergodicity. We note here that, starting from an ansatz very 
similar to the one above, it is possible to obtain the quantum transport equation [14| where 



it is important to relate a given quantum state with the admissible energy surface in phase 
space; thus the above ansatz is in conceptual agreement with the ergodic aspect of many- 
body system. Moreover, this choice fixes the Thomas-Fermi density of states naturally. It 
now becomes important to emphasize that we must restrict ourselves to dilute gas of hard- 
spheres and also assume that the size of sphere is much lesser than the thermal de Broglie 
wavelength. Thus, the ansatz establishes, in fact, a link between RMT and statistical me- 
chanics. We now incorporate the case of two dimensions which otherwise presents enormous 
difficulties. 

In two dimensions, the solutions of the Schrodinger equation, ip(xx, x 2 , • ••,2^)5 under 
an exchange of two coordinates of particles satisfies 

1p(x\, ■ ■ ■ , Xi, ■ ■ ■ , Xj, ■ • ■ , Xn) 

= e i7 "V(xi, • • • , Xj, ■ ■ ■ , Xi, ■ ■ ■ , x N ) (3) 

where v is arbitrary and defines statistics. For v = and v = 1, with (0), one gets the Bose- 
Einstein and Fermi-Dirac distributions. This non-trivial phase and the resulting boundary 
condition arises from the fact that the effective configuration space, Mj^ has a fundamental 
group, 7Ti(M^) = E>n U0|, the Braid group of iV objects which is an infinite, non-abelian 
group. B N is generated by (N — 1) elementary moves o\, • • • , oat-i satisfying the Artin 
relations, 

OiOi+xOi = (Xj + i(XjCXj + i (i = 1, 2, • • • , N - 2) 

OjOi = OiOi, \i- j\ > 2 (4) 



the inverse of cxj is a^ 1 , the identity is denoted by I, and the centre of B n is generated by 
(o" 1 cr 2 • • ' Ojv-i) ■ The multivaluedness of the eigenfunction originates from the phase change 
in effecting an interchange between two coordinates x\ and xf^ (superscripts refering to 
components) which can be expressed as 



V = exp(ii/ i 



i<j 



&j = tan_1 (-7i) fry)- (5) 

The description adopted by us here is referred to as the Anyon Gauge. It is important 
to realise that a set of coordinate configuration can be reached starting from some initial 
coordinates of X particles in an infinite ways, each possibility manifested by an action of an 
element (3 G B N . 

The connection between initial and final sequences is given by ([3]), via the character 
x(P) °f the specific element. Thus, to every (3 G B N , we can associate the affected partial 
amplitude ifj a (P '■ %) jl6| - With one-dimensional unitary representation of the braid group, 



the rudiments of quantum mechanics allow us to write 

$ a (X) = J2 x(P)MP ■ *) (6) 
PeB n 

where ip a (P : X) is the probability amplitude associated in changing a configuration X to 
(P : X) - a configuartion after the action of (3 on X. The wavefunction <& a (X) is to be 
understood as appropriately normalised. The ansatz for Vi[) a ((3 : X) is now 

Vtfj a {(3 :X)=N a J d 2N PA a {(3 : P)5(P 2 - 2mE a )e^' p (7) 

— * 

with A a {(3 : P) satisfying 

(K{h : ft : A)) ME = S j 2 "^ 1 :,^-^^^ , (8) 
(Pi, P2 ^ B N ), and A a (P) satisfy the twisted boundary conditions, 



A a (pi, - ■ ■ ■ ■ ,Pj, - ■ ■ ,p N ) 
e tnu A a (pi , ■ • • , pj, ■ ■ ■ , p u ■ ■ ■ , p N ) (9) 



The question now is in specifying exactly what the matrix ensemble is in this case? The 
form of (|D with A^s not restricted to real, takes into account the T-breaking, and (Q) makes 
the ensemble handed or chiral as a result of P-breaking. Thus (0)-® gives the complete 
description and the ME is, in fact, the chiral- Gaussian Unitary Ensemble (ch-GUE) |17|]. It 
can be easily shown that the Wigner distribution is 

- 1 mL 2 E a N 
Ha ~ N\T(N)E a [ 27rh 2 ] 1 ] 

For the momentum distribution, we need to evaluate the ME-average of $* (P)$*(P") 
with $ = V§ a . With the above ansatz and conditions supplementing it, this average is 



^(p) = fc(p)$;(p")> 



ME 



h 2N 5 a ,,N a ,N, £ £ J2 X*(Pi)x((32)S(P 2 -2mE a 

n,m=0 /3i(m) p2{n) 

(m n ^ 

n - ikv, , .. 



(11) 



where 



Sl N (Q) = h' 2N d 2N Xexp(^Q-X); (12) 

«/ Domain, T> fl 

— * — * 

P is identified with P" after the sum is performed. 
With (|TTD , the momentum distribution is given by 

m) = fdp 2 ...dp N HP) (13) 

/ dpi • • • dp N F(P) 

which formally completes the deduction. However, an exact evaluation of this is very difficult 
and the difficulty is coming from counting of irreducible words formed by the a's. To make 
the precise connection, we give results upto 0{% 2 /L 2 ), an order that is enough for second 
virial coefficient. Leaving the tedium of details to a later publication, we just give our result 
incorporating all elements leading to 0(h 2 /L 2 ), 

F 2 (p t ) = (27rmkT)-> exp ( - 



where 



, !/V , gg=o Sr = + - 2m - 1(e , en ) QP (N) cosjnKu) 
&(N,is) = ^ , (15) 

1 + 2 £~ =1 £^_ 2m(odd) 4 m) (iV) cos(Tri^) 



Qx is the number of elements in Bn composed of 'ml generators whereby the momentum p\ 
is interchanged with another momentum yielding a character exp (i7r.fr z/) (or exp(— iirKv) 
since Qk(N) = QLk(N)); P^ 1 \n) is the number of elements in contributing to 
identity with a character exp(«7r.KV) ( or exp(— iirKv)). Temperature is introduced above 
via the ideal gas law, E a = NkT a . Unfortunately though, this counting problem stands open 
today \TE\. It is very important to note that the ansatz (|?P-([5J) for the special case when 



of = 1 for all i where B^ reduces to symmetric group, Sn, the well-known Fermi-Dirac and 
Bose-Einstein distributions follow. On evaluating pressure, II from (14), denoting area of 
the enclosure by A, we get UA/kT = 1 - (2A)~ 1 X 2 G(N, v), with A 2 = h 2 ^nmkT)- 1 . We 
immediately see that G(N, 0)/(2N) and G(N, 1)/(2N) are 2~ 3 / 2 and -2" 3 / 2 respectively 
yielding the second virial coeffiecient for the Bose and Fermi gases |19[]. For the fractional 
case, with v = even number, 2j + 5 ("boson-based anyons"), comparing our result with ||, 
we get the Sum Rule mentioned in the introduction: 

- 2- 3/2 N- 1 G(N, u)X 2 = (-1 + A\5\ - 25 2 )A 2 /4, (16) 

the right hand side belongs to ||. It is important to note that our deduction is non- 
perturbative and in principle, we can get expressions for higher-order virial coefficients also 
|20| . To understand this, we observe that the relation (8) connects two momentum configura- 



tions of N particles, and not just the momenta of two particles. Thus, it contains information 
that can lead to all virial coefficients. For example, for the third virial coefficient, we need 
to evaluate contributions to F(pi) when three momenta out of N are interchanged. The 
denominator of (13) contains those interchanges which braid three strands in such a way 
that the initial configuration of momenta is preserved whereas the numerator of (13) con- 
tains those which exchange the momentum assignment on all three strands. We have done 
the calculation and the third virial coefficient is expressible in terms of the specific counting 



problem of B^. Here, in order to convince the reader, it suffices to make a comparative 
discussion with the existing calculation. For this, we write down the total contribution to 
the momentum distribution due to a triple interchange emerging from the elements of Bn 
formed by M generators, 

/ fx 2(JV-2) 2A/ 9 

U) E ^l2(N-3){2mE a -^)co6{irku)n^{N) 

\nj _ 2M 6 

+ 2I 2(N _ 3) (2mE a - 3p?) cos(7rA;*/)«Sf (N), (17) 

where Sjf(N) (lZjf(N)) are the number of elements of B^ that (do not) change the mo- 
mentum pi. Id{x) denotes the volume of a D-dimensional hyper-sphere of radius x. The 
reason we give this result here is to show that (17) is a Fourier series with harmonic terms 



like cos2-7tz/, cos4-7rz/, and so on, in complete agreement with the conjectured form |21 |. It is 



becoming evident from the Monte Carlo calculations |22| that the third virial coefficient is 
a series with terms as sin 2 7rz/, sin 4 7rz/, and so on. Our formal result is thus in consonance 
with these works. Also, we mention that (16), (17) and the Monte Carlo estimates provide 
a non-trivial hint on the counting problem itself. 

To summarise, we have formulated in this Letter a way to approach Statistical Mechanics 
from RMT through the ansatz for eigenfunctions which is an essential dynamical input. 
Moreover, we have presented a first-principles-evaluation of the momentum distribution 
of a 'Hard-Sphere Anyon Gas' in the low-density regime which is of a great interest in 
current literature. The nature of difficulties, albeit well-known, are made explicit here in 
direct connection with the existing wisdom (or ignorance) of theory of braid groups. The 
parameter, v has an analogous partner in quantum chromodynamics || and we conjecture 
that the anyon gas discussed here and the z/7r-parametrised quantum chromodynamics belong 
to the same universality class of chiral- Gaussian Unitary Ensemble of RMT. 



S 



I. ACKNOWLEDGEMENTS 



Authors express their deep debt of gratitude to Pierre Gaspard for stimulating discussions 
and constructive criticism. S.R.J, is financially supported by the "Communaute Francaise 
de Belgique" under contract no. ARC-93/98-166. 



9 



REFERENCES 



[1] O. Bohigas, M.-J. Giannoni and C. Schmit, Phys. Rev. Lett., 52, 1(1984). S.W. Mc- 
Donald, Ph.D. Thesis (UMI 8413506, Ann Arbor, MI, 1983). 

[2] S.W. McDonald and A.N. Kaufmann, Phys. Rev. A 37, 3067 (1988); D. Biswas and 
S. R. Jain, Phys. Rev. A 42, 3170 (1990); V.N. Prigodin. N. Taniguchi, A. Kudrolli, 
V. Kidambi, S. Sridhar, "Spatial Correlation in Quantum Chaotic Systems with Time- 
Reversal Symmetry: Theory and Experiment", preprint (1995); M. Srednicki, Comment 
on "Spatial.." by Progodin et al. preprint (1995). 

[3] M. V. Berry, J. Phys. A10, 2083(1977). 

[4] M. Srednicki, Phys. Rev. E 50, 888 (1994). 

[5] A. Lerda, 'Anyons', Lecture Notes in Physics (m-14) (Springer, N. Y., 1992). 

[6] R.E. Prange and S. M. Girvin, Eds. "The Quantum Hall Effect" (Springer, N. Y., 1990); 
R. Ladbury, Phys. Today, July 1995, p. 19. 

[7] J. M. Leinaas and J. Myrheim, Nuovo Cim. 37B, 1(1977); F. Wilczek, Phys. Rev. Lett. 
48, 114(1982); 49, 957(1982). 

[8] D. Arovas, R. Schrieffer, F. Wilczek and A. Zee, Nucl. Phys. B251, 117(1985). 

[9] E. C. G. Stueckelberg, Helv. Phys. Acta 32, 254(1959); 33, 727(1960). 

[10] F. J. Dyson, J. Math. Phys. 3, 1199(1962). 

[11] J. Frohlich and P. Marchetti, Comm. Math. Phys. 121, 177(1989). 
[12] G. Frobenius, J. reine u. angew. Math. 84, 59 (1878). 

[13] C. E. Porter, Statistical Theories of Spectra : Fluctuations, (Academic Press, New York, 
1965), p. 64ff. 

[14] L. van Hove, Physica 21, 517(1955); 25, 268(1959). 

10 



[15] E. Artin, Ann. Math. 48, 101(1947); F. Bohnenblust, Ann. Math. 48, 127(1947). 

[16] M. G. G. Laidlaw and C. DeWitt-Morette, Phys. Rev. D3, 1375(1971); Y.-S. Wu, Phys. 
Rev. Lett. 52, 2103(1984). 

[17] This terminology is used in recent literature, see e.g. J. Verbaarschot, Phys. Rev. Lett. 
72, 2531(1994); we have not examined a relation between our definition and the sense 
in which it is used in this work. We are grateful to Pierre Gaspard for bringing this 
work to our attention. 

[18] J. S. Birman and A. Libgoher, Eds. Contemporary Mathematics : Braids, Proc. Summer 
Research Conference, July 13-26, 1986, Volume 78 (American Mathematical Society), 
see particularly pp. 573-574. 

[19] K. Huang, Statistical Mechanics (John Wiley and Sons Inc., New York, 1987). 

[20] D. Alonso and S. R. Jain, in preparation. 

[21] J. Myrheim and K. Olaussen, Phys. Lett. B299, 267 (1993). 

[22] S. Mashkevich, J. Myrheim, and K. Olaussen, preprint (1996). 



11 



